arXiv:hep-ph/9706452vl 23 Jun 1997 


Interpretations of J/'ip suppression 


N. Armesto and A. Capclla 

Laboratoire de Physique Theorique et Hautes Energies:) 

Universite de Paris XI, Batiment 211, F-91405 Orsay Cedex, France 

Abstract. We review the two main interpretations of J/ip suppression proposed 
in the literature. The phase transition (or deconfining) scenario assumes that 
below some critical value of the local energy density (or of some other geometrical 
quantity which depends both on the colliding systems and on the centrality of the 
collision), there is only nuclear absorption. Above this critical value the absorptive 
cross-section is taken to be infinite, i.e. no J/ip can survive in this hot region. In 
the hadronic scenario the J/ip dissociates due both to nuclear absorption and to its 
interactions with co-moving hadrons produced in the collision. No discontinuity 
exists in physical observables. We show that an equally good description of the 
present data is possible in either scenario. 


1. Introduction 

Charmonium suppression due to Debye screening in a deconfined medium was 
proposed in 1986 by Matsui and Satz |l| and found experimentally by the NA38 
collaboration |2|. However, it was claimed very soon [BJ [4]. [5| that this phenomenon, 
which is also present in pA collisions, could be due to the absorption of the pre¬ 
resonant cc pair in the colliding nuclei. It is nowadays known 0-0 that the data on 
pA and on AB collisions with a light projectile can indeed be described by nuclear 
absorption with an absorptive cross-section a a b s = 7.3 ± 0.6 mb |J. 

Recently the NA50 collaboration has found an anomalous J/ip suppression 
in PbPb collisions, i.e. a suppression which is substantially stronger than the 
one obtained from nuclear absorption with the above value of a a bs [0 0- Two 
different interpretations of this anomalous suppression have been proposed in the 
literature. One is a conventional hadronic interpretation in which there is an extra 

f Laboratoire associe au Centre National de la Recherche Scientifique - URA D0063. 


J/ij} suppression due to final state interaction of the resonant J/'ijj state with co¬ 
moving hadrons B0, |TT| , [T~2[ | . The other interpretation jfj. [bj|. [14], [L5| assumes that 
when the local hadronic energy density is larger than some critical value (taken to 
be around the one reached in a central SU collision), there is a discontinuity in the 
J/-0 survival probability. (See also (TB| for ideas based on percolation of strings.) 


2. Nuclear absorption 


Nuclear absorption is present in all interpretations. We describe it in the 
probabilistic model of Ref. [|J. Let us consider first proton-nucleus collisions. In 
this model, the pre-resonant cc pair is produced at some point z inside the nucleus 
and scatters with nuclei on its path at z' > z, with an absorptive cross-section a abs . 
This produces a change in the A dependence of the J/ip inclusive cross-section. For 
nucleus-nucleus collisions this change is given || by 


S abs (b,s ) 


[1 - exp(-A T a (s) q-qbs)][l - exp(—£> T B (b - s ) a abs )] 

AB T A (s) T B (b - s)al bs j 


Here T A and T B are the nuclear profile functions, determined from a standard Saxon- 
Woods density, and a abs is the absorptive cross-section. In the following we take 
&abs = 7.3 ± 0.6 mb which gives the best fit to the pA data 0. Note that S abs = 1 
for a^s = 0. Expression (|l|) has the meaning of a survival probability of the J/'ip 
due to nuclear absorption. 

Since we are aiming at a quantitative analysis it should be emphasized that the 
probabilistic formula, with its longitudinal ordering in z, can only be true in the low 
energy limit. Therefore it is important to evaluate the uncertainty resulting from 
using this formula at a/s ~ 20 GeV. In a recent paper |l7j the equivalent of Eq. ([]]) 
has been derived in a field theoretical approach. The obtained formula is valid at all 
energies and coincides exactly with ([!]) in the low energy limit. In the asymptotic 


limit we find 17 


S%JI>,s)=exp 


~2 (tAt a{s) 


exp 


-~crBT B (b - s ) 


( 2 ) 


where a is the cc—N total cross-section. At asymptotic energies, the results obtained 
from Eqs. (jTJ) and (|2|) differ only by 8 %. It is amazing that at a/s = 20 GeV the 





























difference between the result obtained with the exact formula, valid at all energies, 
and the one obtained from Eq. (!!]) is less than 1 %. 


3. Phase transition scenario 

The anomalous J/ip suppression observed by the NA50 collaboration in PbPb 
collisions came as a surprise. Indeed, it was believed that the energy densities in SU 
and PbPb collisions were comparable - the larger energies reached in PbPb being 
compensated by a larger interaction transverse area. However, it was soon realized 
that the local energy density (i.e. the energy per unit of transverse area d 2 s ) is larger 
by about 30 % in central PbPb than in central SU collisions [IT3|. More precisely, 


let us consider the well known geometrical factors [|18l [IS 


m A (B)(b, s ) = A(B) T A (b){s) 1 - exp (- a pp B(A) T B[A) {b - s)) . (3) 


The average number of participants is obtained as 

1 


n A + n B = 


&AB 


d 2 b J d 2 s N w (b, s), 


( 4 ) 


where 


N w (b, s ) = m A (b, s ) + m B (b, b- s). 


( 5 ) 


In Ref. |H|, the quantity in Eq. (|5|) is taken as a measure of the energy per unit 
transverse area at each impact parameter. It is found |L3| that the maximum value 
reached in SU collisions at b = 0 is n c = 3.3 fm -2 (see Fig. 1). The model of 
Ref. 0 can then be formulated as follows: For N w (b, s ) < n Cl one assumes nuclear 
absorption alone with a standard value of a a bs (see Section 2). For N w (b , s) > n c , one 
puts cr^s = cxd, i.e. one assumes that none of the J/-0’s produced in this hot region 
can survive. This model gives the maximal suppression that can be obtained in PbPb 
collisions when one imposes that, up to central SU collisions, there is only nuclear 
absorption. The interesting result is that with this simple model one approximately 
obtains the suppression observed experimentally in central PbPb collisions. In Ref. 
[BJ a more detailed analysis of the data has been performed in a similar framework. 
























More precisely, instead of the local energy density N w (b,s), the authors of Ref. || 
consider the quantity 


n(b, s ) 


N c (b,s) 
N w (b,s ) 


( 6 ) 


where 


N c (b, s) = AB T a (s) T B (b — s)a pp . (7) 

The quantity ([61) gives the average number of collisions per participant at fixed b and 
s. It increases with the centrality of the collision. The value of n(b, s) determines 
the onset of deconfinement in Ref. |J . More precisely, at a given impact parameter 
b, the absorptive cross-section is assumed to be infinite for n(b, s ) larger than some 
critical value k c , while for n(b, s ) < n c it takes the standard value a a bs = 7.3 ±0.6 
mb discussed in Section 2. Actually in Ref. [|| one has two distinct values of n c . 
One is for y’s and is taken to be equal to the maximum value of n(b, s) in a central 
SU collision, i.e. 


= Mb = 0 ,s = 0)] sr/ ~ 2.3 . 


( 8 ) 


The second one is for direct J/ij) and is taken as a free parameter satisfying 
nf > k*. The overall J/ip survival probability is then obtained by combining 40 % 
X suppression with 60 % suppression of direct J/^. 

In this way a reasonable description oi J/i\) suppression is obtained. The results 
for PbPb are shown in Fig. 2. 

Finally in Ref. |14| one combines nuclear absorption, deconfinement and 
absorption by co-movers (see Section 4). The quantity which determines the onset of 
deconfinement depends not only on b and s but also on the longitudinal coordinate 
z. When some critical value of this quantity is reached the involved cross-sections 
present a discontinuity but their values remain finite. 




4. Absorption by co-moving hadrons 


The survival probability of the J/-0 due to absorption with co-moving hadrons 
is given by (see [9. (hj, |1^] and references therein) 

'Ny°(b, s) 


S co (b, s ) = exp 


-<j co N™(b,s) In 


N f 


0(N?(b,s)-N f ) 


.(9) 


Here Ny°(b,s) is the density of hadrons per unit transverse area d 2 s and per unit 
rapidity at impact parameter b. Nf is the density at freeze-out that we take to 
be universal. The argument of the log is the interaction time of the J/-0 with co¬ 
moving hadrons. The 0-function is numerically irrelevant (see Section 5). cr co is the 
co-mover cross-section properly averaged over the momenta of the colliding particles 
(the relative velocity of the latter is included in its definition). We treat cr co as a 
free parameter |j20| , |2l|j . All species of hadrons are included in Ny°. This quantity 
has been computed in the dual parton model (DPM) [|22|. It is expressed as a linear 
combination of the geometrical quantities defined in Eqs. (H) and (0) 

Note that S co (b,s ) = 1 for a co = 0. The effects of the co-movers in proton- 
nucleus collisions turn out to be negligeably small (see Section 5). 

The inclusive cross-section for J/ip production in nuclear collisions is then given 
by 



itsib) = — [d 2 s N c (b, s ) S abs (b, s ) S c °(b , s), (10) 

(Jpp J 

where N c is given by i). We see that for Drcll-Yan pair production (a abs = a co = 0), 
1 4 £ = I» y n AB. We take a pp = 30 mb. 

The results Jl^| are presented for three sets of parameters. Set I: a abs = 7.3 mb, 
a co = 0; Set II: cr abs = 6.7 mb, a co = 0.6 mb, Nf = 1.15 firm 2 ; and Set III: o abs = 
7.3 mb, a co = 1.0 mb, Nf = 2.5 fnm 2 . 

The results for J/-0 suppression versus AB are presented in Fig. 3. In Fig. 4 
we show the ratio J/ip over DY versus Et and in Fig. 5 we show the same ratio 
versus L. This variable, defined in Refs. [Q. [8|, is a measure of the centrality of the 
collision. We see that nuclear absorption alone (Set I) fails very badly, whereas Sets 
II and III give a reasonable description of the data. 








5. Proton-nucleus collisions 


Let us first discuss the role of the 9 -function in Eq. (|J). This function looks like 
a discontinuity in the co-mover contribution. However, numerically it has practically 
no effect in our results for SU and PbPb collisions - especially with the parameters 
of Set II. Indeed, the value Nf = 1.15 firm 2 in this Set (which is the same used in 
Ref. 0) is just the average density of co-movers in pp (given by [l/(irRp)]dN/dy). 
Therefore the argument of the log is always larger than one (except for values of the 
integration variables where the argument of the exponent is very close to zero) and 
the 9 -function is irrelevant. With Eq. (|9j), the effect of co-movers is rather small in 
SU (see Figs. 3 and 4) and negligeable in pA. 

In Ref. [|TI1| a different expression for the argument of the log was introduced. 
Using it in Eq.(p|) for AB collisions our results are practically unchanged. However, 
the effect of the co-movers in pA : although small, is then non-negligeable. When 
discussing the ip/DY ratio this effect can be compensated by a decrease of a a b s . 
However, an effect is left in the ratio 0 '/ 0 , which will decrease by roughly 10 % 
between pp and pPb. 


6. 0' suppression 

There is a rather general consensus that (most of) the observed 0 ' suppression 
in SU and PbPb collisions is due to co-mover interactions. Even the authors 0, 
who strongly support the deconfinement scenario for J/0 suppression, do interpret 
in that way the 0 ' over Drell-Yan ratio in all systems (including PbPb). Let me 
discuss this ratio within Set II. Here the value of Nf is the same used in Ref. |6|] and 
the analysis therein is essentially unchanged in our approach. More precisely, in pA 
collisions this ratio is constant both in Ref. [p] and in our formalism (Eq. (|9j)). In 
SU collisions one can choose the value of > af Q such as to reproduce || |10|] the 
ratio 0'/ DY. In the formalism of Ref. [;P] one then gets a good description of PbPb 
using this value of In our case, the value for the last bin of PbPb is too small. 
However, this discrepancy is at the level of two standard deviations, and, therefore, 



one has to wait for better dataf. 

Let me discuss now the effect of deconfinement on ip' suppression. In Ref. 
Tj it is shown that deconfinement produces a sharp decrease of the ratio ip'/ip 


at the critical value of the ip' deconfining phase transitions (taken to be the same 
as for x’s, Eq. (15])). At the time Ref. El was published, it was assumed that 


this deconfinement would take place after pPb. Thus, the decrease of the ratio 
ip'/ip observed in SU was regarded as a qualitative success of the phase transition 
scenario. Nowadays, we know that this phase transition can only take place after 
the last E t bin of SU (Eq. (|8|)). No sharp decrease of the ip'/ip ratio is observed 
here. Therefore the experimental behaviour of the ratio ip'/ip (or ip'/DY) does not 
confirm the qualitative expectations of the phase transition scenario. 


7. The inverse kinematic experiment 

Perturbative QCD calculations yield a very small value of the co-mover cross- 
section, i.e. of the cross-section of a bound J/ip with hadrons near threshold p1[| . 
When non-perturbative effects are introduced, a much larger value is obtained [fn| . 
In order to settle this important question Kharzeev and Satz have proposed to 
measure the J/ip suppression in pA collisions in the backward hemisphere (or A on 
p in the forward one). When the J/ip is slow in the rest system of the nucleus it will 
be produced inside the nucleus and if its interaction cross-section is very small there 
will be no J/ip suppression. The problem, however, is that it may be experimentally 
very difficult to distinguish a cross-section of, say, less than 0.1 mb ]2(J from the 
moderate values (0.5 to 1.5 mb) needed in the co-mover scenario. 

Perhaps a more interesting possibility would be to consider SU in the backward 
hemisphere (or U on S in the forward one). Indeed, in the co-mover scenario the 
decrease of the absorptive cross-section discussed above may be over-compensated 
by an increase of absorption due to an increase of the density of co-movers when one 
moves in rapidity towards the maximum of dN/dy. The observation of a maximum 
in the J/ip suppression at this value of y would support the co-mover scenario. 

f In Ref. fie]] this discrepancy was taken seriously and solved by introducing the exchange reactions 
ip + 7r — > ip' + X and ip' + tt —> ip + X. 






A further advantage of this proposal is that it can be done in the rapidity range 
(0 < y C m < 1) covered by the present NA50 experimental set np. The calculation 
of this effect is, of course, possible in the co-mover scenario but it requires a 
hypothesis on the energy dependence of a a b s . It should be stressed, however, that 
the observation of such a maximum does not rule out the phase transition scenario. 
Indeed, the J/i[) environment at — 1 < y cm < 0 in SU (or at 0 < y cm < 1 in US) 
is hotter than the one in SU at 0 < y cm < 1, and therefore deconfinement can be 
present there. This shows that the geometrical criteria for deconfinement used in 
can only be appropriate for a flat rapidity density of hadrons. When 
this is not the case the critical value of the local energy density has to be a function 
of y. It also shows the difficulty to disentangle the co-mover scenario from the 
deconffiiing one. 

8. Discussion and conclusions 

The main difference between color deconfinement and co-mover absorption is the 
presence in the former of a sudden change in the J/-0 survival probability, whereas 
all changes are smooth in the latter. However, this discontinuity is assumed to take 
place in a variable (the local energy density) which is not directly observable. 

As for the number of free parameters, in the deconfinement approach of Ref. 
PH there is (apart from a a b s which is present in all approaches) only the value of 
n c (the critical value of the local energy density). Moreover, one has chosen a value 
of Cabs = oo above the critical density. In Ref. || there is an extra parameter since 
one takes two distinct critical values for %’s and for direct J/^j deconfinement. 

In the co-mover scenario there are two free parameters: cr co and Nf. The Et (or 
L) dependence in PbPb is controlled by both quantities. However, we see in Figs. 
3, 4 and 5 that changing Nf by a factor 2 has very little effect on the results, since 
this change is compensated by a corresponding change in the value of a co . Because 
the effects of these two parameters are strongly correlated, it is not possible in the 
simple co-mover approach of Refs. 0 . 0.0 to reproduce the structure in the Et 
(or L) dependence of the PbPb data. Taking a mixture of %’s and direct J/tf} with 
distinct parameters, as done in Ref. [6{], would improve the situation. However, one 
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should make sure that such a structure is not due to a statistical fluctuation. 

The success of the co-mover model is due to the fact that, when a couple (Nf, a co ) 
is chosen in such a way to reproduce the Et (or L ) dependence in PbPb , the 
corresponding effect of the co-movers in SU turns out to be rather small and does 
not spoil the success of the nuclear absorption model. Moreover, the dependence on 
centrality of the local density of co-movers is stronger in PbPb than in SU collisions. 
This produces the change in the L slope between the two systems, seen in Fig. 5. 

In conclusion, to distinguish the deconfinement scenario involving a phase 
transition (or discontinuity) from the co-mover scenario turns out to be a 
quantitative rather than a qualitative issue. The present data can be reasonably 
well described in both approaches with a comparable (and small) number of free 
parameters. So far no experiment has been suggested which would allow to 
disentangle these two mechanisms. 
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Figure 1. The density N w (b, s ) in Eq. (||) for s along the direction of the impact 
parameter, for various values of the impact parameter b = 0, 2,4,... fm. The origin 
is at a distance 6/(1 + Rb/Ra) from the center of nucleus A. Left: SU collision; 
right: PbPb collision. The horizontal dashed line coresponds to the largest density 
achieved in the SU system, n c = 3.3 fm^ 2 . 








Figure 2. The experimental J/ip survival probability over 
by the deconfinement suppression in PbPb collisions, with n 



DY ratio [[/) divided 
* = 2.3 and = 2.9 



Figure 3. J/^> suppression versus AB: Set I (dotted line), Set II (solid line) 
and Set III (dashed line) compared to the experimental data (8). Note that the 
calculations have been performed only for those nuclei where data exist. The 
obtained values have been joined by straight lines. For clarity of the figure, the 
results of Set III are shown only for SU and PbPb. 







Figure 4. J/ip over DY ratio versus Et : Set I (dotted line), Set II (solid line) 
and Set III (dashed line) compared to the experimental data Q. 



Figure 5. Same as in Fig. 4 plotted versus L j|, ||]. 







